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ABSTRACT 

We present an analysis of the _B-band and V-band rise-time distributions of nearby 
Type la supernovae (SNe la). Drawing mostly from the recently published Lick Ob- 
servatory Supernova Search sample of SNe la ([Ganeshalingam et al. 20101, together 
with other published nearby SNe la with data starting at least one week before max- 
imum light, we use a two-stretch template-fitting method to measure the rise and 
decline of BV light curves. Our analysis of 61 SNe with high-quality light curves in- 
dicates that the longer the time between explosion and maximum light (i.e., the rise 
time), the slower the decline of the light curve after maximum. However, SNe with 
slower post-maximum decline rates have a faster rise than would be expected from a 
single-parameter family of light curves, indicating that SN la light curves are not a 
single-parameter family of varying widths. Comparison of the i?-band rise-time distri- 
bution for spectroscopically normal SNe la to those exhibiting high-velocity spectral 
features indicates that high- velocity (HV) SNe la have shorter i?-band rise times com- 
pared to their spectroscopically normal counterparts. After normalising the i?-band 
light curves to Ami5(-B) = 1.1 mag (i.e., correcting the post-maximum decline to 
have the same shape as our template), we find that spectroscopically normal SNe la 
have a rise time of 18.03 ± 0.24 d, while HV SNe have a faster i?-band rise time of 
16.63 ±0.29 d. Despite differences in the B band, we find that HV and normal SNe la 
have similar rise times in the V band. Furthermore, the initial rise of a SN la i?-band 
light curve follows a power law with index 2. 20^o'ig, consistent with a parabolic rise 
in flux predicted by an expanding fireball toy model. We compare our early-time B- 
band data to models for the predicted signature of companion interaction arising from 
the single-degenerate progenitor scenario. There is a substantial degree of degeneracy 
between the adopted power-law index of the SN light-curve template, the rise time, 
and the amount of shock emission required to match the data. 
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1 INTRODUCTION 

Type la supernovae (SNe la) are believed to be the result 
of a thermonuclear runaway explosion of a C/O white dwarf 



(WD) approaching the Chandrasekhar limit (see Hillebrandt 
|fe Niemeyer||2000| for a review). Explosive nucleosynthesis 
up to bb Ni releases ~ 10 51 erg, unbinding the progenitor. 
The subsequent light curve is powered by injection of en- 
ergy from the radioactive decay of 56 Ni; the 7 rays degrade 
to longer wavelengths as they diffuse out through the ex- 
panding ejecta. 

The well-established relationship between the light- 
curve width and the luminosity at peak brightness allows 
SNe la to be "standardizable candles" at optical wavelengths 
( Phillips| 1993 1 , and possibly almost standard candles in the 
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infrared (Krisciunas et al. 2004 Wood-Vasey et al. 2008 



Folatelli et al. 20101. Application of the width-luminosity 



relation along with colour information to SNe la over cos- 
mological scales has led to the discovery of the accelerat- 
ing expansion of the Universe ( Riess et al.|1998| |Perlmutter 



2009b Amanullah et al. 



et al.||1999] see also |Astier et al.||2006l |Riess et al.||2007| 
Wood-Vasey et al.|2007||Kowalski et al.||2008| |Hicken et aL] 



2010 1, indicating either the pres- 



ence of "dark energy" having a negative pressure or a failure 
of general relativity on the largest scales. Recent work has 
also shown that including spectral flux ratios may also aid 



in reducing Hubble-diagram residuals (Bailey et al. 2009 
Blondin et"aL]|2011 1. 

Despite the successful cosmological application of 
SNe la, the SN community lacks a clear understanding of the 
nature of their progenitor systems. Possible scenarios include 
a single-degenerate WD paired with a red-giant post-main- 
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sequence star undergoing stable mass transfer until Mwd 
approaches M C h ~ 1.4 M Q dWhelan & Iben||1973| iLivio 



2000 1 , a double-degenerate WD merger that reaches or sur- 
passes M Ch ( |Webbink|[T984l |Iben fe Tutukov||l994] |, or the 
result of a sub-Chandrasekhar explosion of a WD steadily 



burning helium accreted from a companion ( Shen & Bildsten 



2009). Each scenario carries with it tentative observational 



evidence and hardships imposed by unrealised theoretical 
predictions. 

While considerable attention has been focused on the 
post-maximum decline of the SN la light curve, less has 
been paid to the rise of the SN la from explosion to max- 
imum brightness. This is due to the dearth of data in the 
days following the SN explosion and the intrinsic difficulty 
of finding SNe la shortly after explosion. The Lick Observa- 
tory Supernova Search (LOSS) has been successful at find- 
ing SNe in the nearby Universe (redshift z < 0.05) and the 



Sloan Digital Sky Survey (SDSS) Supernova Survey (Frie- 
man et al. 2008) has found and monitored 391 SNe at mod- 
erate redshift (0.03 < z < 0.35). Higher redshift searches 



such as the SuperNova Legacy Survey ( Astier et al. 2006 1 



and ESSENCE ( |Wood-Vasey et al.||2007fr benefit from cos- 



mological time dilation, allowing the SN rise to be spread 
out over more days in the observer's frame. 

A first attempt at quantifying SN la rise times by 
|Pskovskii|(1984) analysed 54 historical light curves gathered 



from photographic plates and visual B magnitude estimates, 
finding a typical rise of 19-20 d. Recent attempts have made 
use of more reliable data taken with charge-coupled devices 
(CCDs) which are linear in translating photons to counts 
over a large dynamic range. A seminal paper by Ries s et al.| 
( 1999b hereafter R99) used observations of nearby SNe la 



to find a B-band rise time of 19.5 ±0.2 d for a normal SN la 
having a decline of 1.1 mag between maximum light and 15 d 
after maximum in the B band (i.e., Amis(_B) = 1.1 mag; 
Phillips 1993). Studies using data on higher redshift SNe la 
have found a rise time in concordance with that of R99 
( Aldering et al.|2000 Conley et al.|2006 1 , advancing the no- 
tion that there is limited, if any, evolution in SN la properties 
from high to low redshifts. Most recently, ( [Conley et al. 2006 1 
hereafter C06) found a rise of 19.10Jlg;^ (stat) ± 0.2(sys) d 
using SNLS data. 

In an analysis of B- and U-band photometry of eight 
nearby SNe la with excellent early-time data, |Strovink| 
( 2007 hereafter S07) found tentative evidence for two popu- 



lations of B-band rise times. After correction for light-curve 
decline rate, S07 found that three SNe rise in 18.81 ± 0.36 d 
and five SNe rise in 16.64 ± 0.21 d. More recently, Hayden 
et al. (2010a hereafter H10) analysed B- and l/-band pho- 



tometry 105 SNe la from the SDSS sample and found that 
SNe la come from a rather broad distribution of B-band 
rise times, with evidence indicating that slower declining 
events (e.g., more luminous SNe) have some of the fastest 
rise times. This result has significant implications for light- 
curve fitting techniques that rely on a single-parameter fam- 
ily of light curves and theoretical modeling of SNe la. H10 
find an average B-band rise time of 17.38 ± 0.17 d, a depar- 
ture from the results of R99 and C06 which H10 trace back 
to differences in fitting methods. 

In this paper, we analyse available data on nearby 
SNe la to measure the rise times, relying heavily on the re- 
cently released LOSS sample (Ganeshalingam et al.||2010). 



Previous analyses such as those of S07 and H10 use combined 
results from both B- and V-band photometry to measure the 
B-band rise time. In this paper, we compare the B-band rise 
time measured in the B and V bands and find that such a 
combination may not be appropriate. Instead, we present 
independent analysis of the B and V bands. Specifically, we 
define the rise time in a photometric band as the elapsed 
time between explosion and maximum light for that partic- 
ular band. Nearby SNe offer the advantage of being able to 
be monitored by small-aperture telescopes and benefit from 
not requiring significant i^-corrections, which at early times 
are ill defined because of a lack of available spectra to model 
the SN spectral energy distribution (SED). 



2 DATA 



2.1 LOSS Light Curves 



LOSS is a transient survey utilizing the 0.76-m Katzman 
Automatic Imaging Telescope (KAIT) at Lick Observatory 



(Li et al. 2000; Filippenko et al. 2001 see also Filippenko, Li, 



& Treffers 2011, in prep.). KAIT is a robotic telescope that 
monitors a sample of ~ 15, 000 galaxies in the nearby Uni- 
verse (z < 0.05) with the goal of finding transients within 
days of explosion. Fields are imaged every 3-10 d and com- 
pared automatically to archived template images, and po- 
tential new transients are flagged. These are subsequently 
examined by human image checkers, and the best candi- 
dates are reobserved the following night. Candidates that 
are present on two consecutive nights are reported to the 
community using the International Astronomical Union Cir- 
culars (IAUCs) and the Central Bureau of Electronic Tele- 
grams (CBETs). Time allotted to our group on the Lick 
Observatory 3-m Shane telescope with the Kast double spec- 
trograph ( Miller fe Stone| | 1993[ ) is used to spectroscopically 
identify and study candidates. Between first light in 1997 
and 2010 September 30 UT, LOSS found over 865 SNe, 382 
of which have been spectroscopically classified as SNe la. 
The statistical power of the LOSS SNe is well demonstrated 
by the series of papers deriving the nearby SN rates ( |Leaman| 
et al.|2011||Li et al.|2011a|b| >. 

In addition to the SN search, KAIT monitors active 
SNe of all types in broad-band BVRI filters. The first data 
release of BVRI light curves for 165 SNe la along with de- 
tails about the reduction procedure have been published 



by Ganeshalingam et al. (20101. In summary, point-spread 



function (PSF) fitting photometry is performed on images 
from which the host galaxy has been subtracted using tem- 
plates obtained > 1 yr after explosion. Photometry is trans- 
formed to the Landolt system ( |Lando lt 1983, 1992) using 
averaged colour terms determined over many photometric 
nights. Calibrations for each SN field are obtained on pho- 
tometric nights with an average of 5 calibrations per field. 

The LOSS light curves represent a homogeneous, well- 
sampled set of BVRI light curves. The average cadence is 
3-4 d between observations, with a typical light curve having 
22 epochs. Of the 165 BVRI light curves in the sample, 70 
have data starting at least one week before maximum light. 
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2.2 Light Curves from Other Nearby Samples 

In addition to the LOSS sample, here we include data from 
the following previously published SN la samples: the Calan- 
Tololo sample ( Hamuy et al.|[l996| l, the Center for Astro- 
physics (CfA) Data Releases 1-3 ( |Riess et al. 1999a Jha 
et al.|[2006a Hicken et al.|2009a[), and the Carneg ie Super- 
nova Project (CSP) dataset ( [Contreras et al.|2010 l. In cases 
where there are data from multiple samples, we chose the 
dataset with the best-sampled light curve to avoid intro- 
ducing systematic calibration error. With the exception of 
the CSP sample, all light curves are in the Landolt pho- 
tometric system. The CSP light curves are in the natural 
system of the Swope telescope at Las Campanas Observa- 
tory. Comparing LOSS B-band light curves (in the Landolt 
system) with CSP B-band light curves (in the Swope nat- 
ural system), we find differences of ~ 0.03 mag which do 
not appear to be correlated with SN colour. We adopt 0.03 
mag as the systematic uncertainty of the CSP Swope system 
light curves. 

We also include individual light curves for SN 1990N 
( |Lira et al.|p98]) SN 1 992A ( jAltavilla et~aL] [2004]) SN 
1994D ( jPatat et al.|1996|, SN 1998aq ( |Riess et al.pX)5]), SN 
1999e e ( [Stritzinger et al.|2002||, SN 2 003du ( |Stanishev et al. 
2007| ), SN 2007gi ( |Zhang et al.|2010[ ), and a preliminary re- 
duction of SN 2009ig which was found by LOSS about 15 d 
before maximum light. In total, we have BV light curves for 
398 SNe. 



3 METHODS 

In this section we detail the method used to measure the 
rise times of our sample. Most previous measurements of 
the SN la rise time have used a single-stretch fit (see R99, 
C06). These studies "stretched" the template light curve 
along the time axis to determine the stretch value, s, that 
best fit the data. A light curve narrower than the template 
would have s < 1, a wider light curve s > 1, and a light 
curve that matched the template perfectly s = 1. Explicit in 
this formalism is that a single stretch value applies to both 
the rising and falling portions of the light curve. Instead, 
we adopt a two-stretch fitting procedure first introduced by 
H10 to fit our B- and V-band data; we fit each band inde- 
pendently. Here we discuss our implementation of the two- 
stretch fitting routine with a template created mostly from 
LOSS data. 



3.1 The Two-Stretch Fitting Method 

In our two-stretch fitting routine, pre-maximum and post- 
maximum data are decoupled, allowing the two portions of 
the light curve to take on different stretch values to match a 
template. We define s r to be the "rise stretch" and s/ to be 
the "fall stretch." The template has s r — Sf = 1.00 by con- 
struction. The two stretched light-curve portions are joined 
at peak where the first derivative is 0, ensuring a continu- 
ous function with a continuous first derivative at maximum 
light. 

Mathematically, the time axis is stretched such that 



Sr .(i+ 2 ) 

t-tg 
s f (l+z) 



t > t , 



where r represents the effective "stretch-corrected" rest- 
frame epoch, to is the time of maximum light, and z is the 
SN redshift. 

Similar to other implementations of stretch (e.g., |Gold-| 
haber et al.|2001 l, we model each light curve as a function 
of time by 



/(*) = /oS(r), 



(2) 



where /o is the peak flux and S is our normalised template 
light curve. We perform a x 2 minimization to the quantity 
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where i is an index summed over all observations, F(ti) are 
individual flux measurements, <7 p h t is the photometric un- 
certainty including systematic error, and ^temp late is the un- 



certainty in our template as described in [ 3.2 Our \ mini- 



mization fits for the values of /o, to, s r , and s/. We restrict 
the fit to data within r < +35 d relative to maximum light, 
after which SNe enter the nebular phase where the stretch 



parametrisation is no longer applicable ( Goldhaber et al 
200Tj ). 

The rest-frame rise time of a light curve (i.e., the elapsed 
time between explosion and maximum light in that band) 
is obtained by multiplying the measured s r by the fiducial 
rise time of the template. Following S07, we define the B- 
band fall time to be the amount of time required for the B- 
band light curve to decline by 1.1 mag starting at maximum 
light in B. We define the fall time for the V band to be the 
required time for a l^-band light curve to fall by 0.66 mag; 
the fall time of a light curve is 15 Sf d. By construction, both 
templates with s/ = 1 have a fall time of 15 d. 

The _B-band fall stretch s/ is directly related to 
Ami5(B) by the simple, empirical formula 

Am 15 (B) « l.l-1.70(s/(B)-1.0)+2.30(s/(B)-1.0) 2 , (4) 

which was found by stretching our B-band template and 
reading off the resulting value of Amiz>{B) mag. Similarly, 
for our V-band template, we find 

Am 15 (V) « 0.66 - 0.83(s/(V) - 1.0) + 0.94(s/(y) - 1.0) 2 . 

(5) 
Note that Equations [4] and [5] are only valid for the LOSS 
templates, and that Sf(B) — 1 and Sf(V) — 1 respectively 
correspond to Ami5(B) = 1.1 mag and Ami5(V) = 0.66 
mag. 

Jf-corrections are computed with the spectral series of 



Hsiao et al. ( 2007 1 which provides the spectral evolution of 



(1) 



a SN la with a one-day cadence. Using all available multi- 
colour photometry for a SN, the spectrum for the phase near- 
est the photometry epoch is warped to match the colours in 
the observer's frame using a third-order spline with knots 
placed at the effective wavelength of each available filter, 
^-corrections are computed from the warped spectrum and 
are typically < 0.05 mag. 

For a single SN, our fitting procedure gives the date of 
maximum light, maximum flux, rise stretch, and fall stretch 
for both the B- and V-band light curves. For the purpose 
of comparing our results to those of H10, we compute the 
B-band rise time (i.e., the time between explosion and max- 
imum light in B) found using both B- and V-band pho- 
tometry. To determine the B-band rise time using l/-band 
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photometry, we measure the V-band rise time (the time be- 
tween explosion and maximum light in V) and subtract off 
the time between maximum light in V and maximum light in 
B. We employ this determination of the B-band rise time us- 
ing 1^-band photometry to compare the rise-time behaviour 
of the two bandpasses over the same time period. Conse- 
quently, the determination of the B-band rise time using 
the V band typically has larger uncertainties than using the 
B-band photometry, since we must also include the error 
from both the times of B and V maximum. We emphasise 
that the V-band rise time is defined to be the elapsed time 
between explosion and maximum light in the V band. 



3.2 Template 

The template plays an important role in measuring light- 
curve properties. A template which does not reflect the data 
will lead to fits with systematically incorrect measurements 
of light-curve properties (H10; Aldering et al.| |2000[). Here 
we discuss the sample of objects and the method used to 
construct our light-curve templates for the B and V band. 
We construct our templates from a sample of well- 
observed "normal" SNe la. Most of the objects come from 



the set of SNe la observed by LOSS ( Ganeshalingam et al 



2010), but we also include several objects published previ- 



ously. SNe la with excellent light curves but known peculiar- 
ities such as SN 1999ac ([Phillips et al.|2006| ) SN 2000cx ([Li] 
et al.|2001a| |Candia et al.|2003[), SN 2002cx (|Li et al|2003| 



2007), and SN 2009dc (Yamanaka et al. 2009 Silverman 



Jha et al.||2006bp, SN 2004dt (Leonard et al.||2005| |Wang 
et al.|2006J|Altavilra et al.|2007~SN 2005hk (|Phillips et al 



et al. 20101, are avoided. The underluminous SN 1991bg- 



like objects show distinctly different photometric behaviour 
compared with the rest of the SN la population, so they are 
excluded from the sample as well. In total, the set includes 
60 objects, many of which are also in the rise-time SN la 
sample discussed in |3.4| 

For the zeroth-order template light curve, we adopt 
the fiducial A = (Amis(B) = 1.1 mag) template from 



the Multi-co lour Light-Curve Shape (MLCS2k2) fitter ( |Jha 
et al. 2007 )rj Light-curve data are then fit using a two- 
stretch fitting parametrisation. The parameters being fit via 
a x 2 minimization technique are the time of maximum light 
to, the flux at maximum light /o, the rise stretch s r , and 



the fall stretch Sf (see E 3. 1 for details) after correcting the 
light curves for time dilation using redshifts obtained from 
the NASA/IPAC Extragalactic Database^The light-curve 
data are then normalised using the best-fit parameters to 
have a peak flux of 1.0 at r = and de-stretched, using s r 
and Sf, along the time axis to match the template. After 
normalizing and de-stretching all of our light-curve data to 
match the shape of the template, we study the mean resid- 
uals between the data and the template. 

While we find that even though the A = templates 
do a reasonable job in fitting the data, there are still small 
systematic trends in the fit residuals. For each band, we fit 
a smooth curve to the residuals and use it as a correction 



1 Downloaded from www.physics.rutgers.edu/~saurabh/mlcs2k2/ 
on October 7, 2010. 

2 http://nedwww.ipac.caltech.edu/ 



to the input template. These "refined" templates are then 
used to fit the light-curve data from our template sample 
again and the fit residuals are studied. This process is iter- 
ated until convergence is reached — that is, no systematic 
trend is observed in the fit residuals. Convergence is achieved 
within 5 iterations. We restrict this part of the template 
training procedure to data where the MLCS2k2 template is 
well defined, but before the SN enters the nebular phase. 
For the B band this is within the range — 10 < r < +35 d 
with respect to maximum light in B. For the V band this is 
— 11 < r < +35 d with respect to maximum light in V. 

To estimate the uncertainties of our templates, we bin 
the fit residuals in 3 d intervals and calculate their root- 
mean square (RMS). Because a particular light-curve fit 
can have systematic residuals relative to the input template 
(i.e., several data points in one portion of the light curve 
all show negative residuals, while data in another portion 
all show positive residuals), the residuals at different epochs 
are correlated, and the RMS measurements are an overes- 
timate of the true uncertainties. It is difficult to quantify 
how the residuals are correlated because different light-curve 
fits have different patterns of residuals. We assume that the 
residuals are equally affected by the correlated errors in all 
portions of the template, and apply a constant scaling factor 
to the RMS measurements so that the overall fit to all of the 
data has a reduced \ 2 ~ 1- These scaled RMS measurements 
are adopted as the uncertainties of our templates. 

Suffering from a dearth of early-time data, we adopted 
an expanding fireball model to describe the SN light curve 
for r < — 10 d relative to maximum light based on the ar- 
guments presented by R99 and employed in most rise-time 
studies thereafter. After explosion of the progenitor, the SN 
undergoes free, unimpeded expansion such that its radius, 
R, expands proportionally with time, r. Approximating the 
SN as a blackbody, the optical luminosity through a broad- 
band filter on the Rayleigh- Jeans tail of the SED is given by 
L oc R 2 T oc (vt) 2 T. Despite recent observations indicating 
that T may actually change substantially over this period 
( Pastorello et al]|2007 Hayden et al.pOlOa l, if we assume 
that changes in v and T are modestly small, then L oc r 2 . 
Wrapping our ignorance into a "nuisance parameter" a, we 
can write the flux in the rise-time region as / = a(r + t r ) 2 , 
where t r is the rise time of the template. 

To determine the rise time for our template, we restrict 
our sample to light curves with data starting at t ^ — 10 d 
relative to B-band maximum and at r ^ — 11 d for the V 
band. The following approach was adopted. 



(i) Create a random realization of the light curves in our 
sample using reported photometry uncertainties for each 
data point and systematic calibration uncertainties for each 



light curve in our template section (see [ 3.3 for details on 
how this is implemented). 

(ii) Perform the two-stretch fitting routine, restricting the 
fits to data within the range —10 < r < 35 d for B or 
— 11 <r<35d for V . The fit is restricted to this range to 
avoid imposing a shape in the region where we will fit for 
the rise time. 

(iii) Normalise and two-stretch correct the light curves in 
the sample. 

(iv) Fit a parabola to the ensemble of light-curve points 
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with r ^ — 10 d for B or r ^ — 11 d for V with the constraint 
that the parabola is continuous with the template. 

The procedure outlined above was performed 1000 
times for the B- and V-band data independently. We find 
a best-fit rise time of 17.92 ± 0.19 d for the B band and 
19.12 ± 0.19 d for the V band. 

We also attempted to find the best rise time for our 
template using the approach outlined by H10, except on a 
finer grid. We tested templates with different rise times in 
the range 15 ^ t r ^ 21 d in 0.1 d increments (as opposed to 
the 0.5 d increments of H10). For each template, light curves 
with data starting at least ten days before maximum light 
were fit using the two-stretch method (see E 3. 1 for more de- 
tails). Unlike the previous procedure to determine the rise 
time, all data at r < +35 d are used in the fit. We calculated 
the x 2 statistic for all stretch-corrected data in the region 
r ^ —10 d for each template. We fit a fourth-order poly- 
nomial to the curve of reduced x 2 as a function of t r . The 
minimum of the polynomial is taken to be the rise time that 
best fits our data. 

Performing a Monte Carlo simulation of this procedure 
to draw 500 unique realizations of our dataset, we find a 
best-fit rise time of 17.11 ± 0.09 d for our B-band template 
and 18.07±0.11 d for our V-band template. These rise times 
disagree at the 4-5<r level with what we found above. It is 
unclear why these two procedures find such different results 
for the best-fit rise time for our template. However, when 
comparing a B-band template with a rise time of 17.11 d 
to the data after they have been two-stretch corrected using 
stretches found with fits restricted to — 10 < r±35 d, we see 
a significant systematic trend for data r < — 13 d. Similar 
results are found in a V-band comparison with a template 
rise time of 18.07 ± 0.11 d. While the nature of the discrep- 
ancy eludes us, we use a S-band rise time of 17.92 ±0.19 d 
and a V-band rise time of 19.12 ±0.19 d to avoid intro- 
ducing any systematic error. As a precaution, we have run 
our analysis using both sets of template rise times and find 
that while some of the final numbers change, they change 
systematically (by ~ 0.50-0.75 d) and do not affect any of 
our qualitative results or final conclusions. In §4.5| we dis- 
cuss using different rise times in the context of searching 
for signatures of interaction with the companion star in the 
single-degenerate scenario. 

We caution the reader from interpreting the above tem- 
plate rise times as the "typical" rise time of a SN la or even 
the rise time of a SN with Ami5(-B) = 1.1 mag. The rise 
time found above should be viewed as the rise time required 
to match the shape of our template light curve in the region 
— 10 < t < d. The rising portion (r < d) of our template 
light curve based on the MLCS2k2 template does not a pri- 
ori reflect a light curve with a decline /S.m\^(B) = 1.1 mag. 
Thus far, the goal of constructing our template was to find 
a light-curve shape that will fit our sample by applying in- 
dependent stretches to the rising and falling portions of the 
light curve. To that extent, the actual number associated 
with the template rise time is meaningless; the significance 
is in the shape of the light curve. For example, we could 
construct an equivalent B-band template with a rise time 
of 35.84 d by stretching the rising portion (r < d) of 
our 17.92 d template by a factor of 2. Using a template 
with a rise time of 35.84 d would decrease the measured 
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Figure 1. Stretch-corrected B-band light curves for 61 objects 
using the two-stretch method. All light curves have been shifted 
such that t = d is at maximum light and they have the same 
peak flux. Ovcrplottcd in red is the LOSS template. We restrict 
our fits to data within 35 d of maximum light. The inset panel 
is a detailed view of the early-time rise of the SN. The bottom 
panel shows the residuals between the stretch-corrected data and 
the template normalised by the template flux. 



rise stretches by a factor of 2. The final rise time for each 
SN found by multiplying the rise stretch by the fiducial rise 
time of the template will be the same for both templates. 
In §4.1[ we will use the shapes of our template light curves 
to measure the rise and fall of our data sample and we will 
address the fall stretch-corrected rise time of SNe la in our 
sample. 

Figure [T] compares our template light curve to the 
stretch-corrected light-curve data. Residuals to the fit scaled 
by the template flux are plotted in the bottom panel of the 
figure. The template fits all portions of the light curve, with- 
out any systematic trends in the residuals. 



3.3 Estimating Uncertainties 

To estimate the uncertainties in our fitting procedure, we 
use a Monte Carlo procedure including the effects of system- 
atic calibration error from different photometric surveys, the 
uncertainty in rise time of our template, and the statistical 
photometric error. The prescription for one simulation for 
our Monte Carlo procedure is as follows. 

(i) For each survey, model calibration uncertainties by 
choosing a random photometric offset, and change all pho- 
tometry from that survey by that random amount. The ran- 
dom offset is chosen for each survey assuming a Gaussian 
distribution with a mean offset of 0.0 mag and a = 0.03 
mag. SNe from the same survey will have the same photo- 
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Figure 2. The redshift distribution of our sample broken into the 
three spectroscopic subclassifications: Normal, High Velocity, and 
SN 1991T/1999aa-like. The median redshift for the entire sample 
of 61 objects is 0.016. 



metric offset while SNe from a different survey may have a 
different offset. 

(ii) Model the photometric error by perturbing every pho- 
tometric point for each SN randomly based on its reported 
uncertainty, assuming the uncertainty is Gaussian. 

(iii) Modify the LOSS template to have a rise time given 
by a random draw from a Gaussian distribution defined by 
the mean of 17.92 d and a = 0.19 d for B or 19.12 d and 
a = 0.19 d for V (see j ]3.2| for how these values were derived). 

(iv) Fit simulated data with the modified template. 

We perform 1000 Monte Carlo simulations for our B- and 
V-band photometry independently. We take the parameter 
uncertainties to be the standard deviation in our 1000 trials 
and keep track of the covariance matrix between fit param- 
eters. 



3.4 The Rise-Time Sample 

We restrict our final sample to objects that have at r f ^ 

1.5 d in both the B and V bands to ensure that we are us- 
ing objects where the rise and fall of the light curve are being 
fit well. By the nature of our Monte Carlo procedure, each 
light curve does not have a single value of \ 2 associated with 
it. We instead identify SN light curves which have a median 
reduced \ 2 < 1-5 over the 1000 trials. Fits are visually in- 
spected to identify the poor ones, which are excluded from 
our sample. We also place the somewhat strict requirement 
that the SN have data starting 7 effective days before maxi- 
mum (defined by EquationfTl to anchor the measurement of 
s r . We explore the effects of relaxing or tightening these re- 
quirements in §4.7| and find that changing the requirements 
does not affect our final results. 

Objects similar to the subluminous SN 1991bg ( |Filip-| 
penko et al.||1992b Leibundgut et af||1993 1 were unable to 
be fit satisfactorily by our fitting routine. This is not sur- 
prising given the significant differences in light-curve shape 
and spectral evolution between the SN 1991bg-like SNe la 




10 20 

Days Past Max 



Figure 3. B-band light curves for the SNe used in this analy- 
sis. The light curves have been shifted such that they share the 
same peak flux and t = is the time of maximum light. The light 
curves are well sampled and start at least one week before maxi- 
mum light. The SNe are coded by Amis(B) which measures the 
the post-maximum decline. The collection of light curves strongly 
suggests that SNe with slow rise also have a slow decline. 



and "Branch-normal SNe la" ( Branch & Tammann 1992 1 
In " 



4.6| we discuss our attempt to measure the rise times of 
SN 1991bg-like objects. We also exclude SN 2005hk, which 



2003 Jha 



is of the peculiar SN 2002cx subclass (Li et al. 
|et al.|2006bl|Phillips et al.|2007| >. 

Of the initial BV light curves for 398 SNe, 95 have data 
starting at least 7 effective days before maximum light that 
could be fit with our two-stretch routine. Of the 95, 63 have 
good fits with rise and fall times that be can reasonably 
measured to within an uncertainty of 1.5 d. 

For the purposes of analysing the rise time of differ- 
ent spectroscopic subclasses, we further break our sample 
into three groups: Normal, High Velocity (HV), and SN 
1991T/1999aa-like (|Filippenko et al.|1992a||"Li et al.|2001b 
Garavini et al. 2004). For the classification of HV objects, 



we adopt the criterion of Wang et al. ( 2009 1 that the average 
vsi for spectra taken within one week of maximum light is 
> 1.2 x 10 4 km s _1 . Objects which have usi ^ 10 4 km s _1 
(more typical Si velocities for a normal SN la) are classified 
as normal SNe la. 

We also spectroscopically identify SN 1991T/1999aa- 
like objects in our sample. The combination of SN 1991T-like 
and SN 1999aa-like objects is based off of previous stud- 
ies which note the similar photometric and spectroscopic 
properties of the two subtypes (|Li et aL 2001b Strolger 



et al. 20021. In particular, both exhibit broad light curves 
(Ami5(_B)« 0.8-0.9 mag) and spectra indicative of higher 



photospheric temperatures than normal SNe la (see Filip- 
penko 1997 for a review of the spectroscopic diversity of 
SNe la). 

Our subclass information comes from spectra taken 
within one week of maximum light observed as part of the 
Berkeley SuperNova la Program (BSNIP; Silverman et al. 
2011, in prep.). Using a modified version of the SuperNova 
IDentification (SNID) code ( |Blondin fc Tonry|2007J ), we are 
able to classify SN spectra using a cross-correlation algo- 
rithm against a spectral database of know subtypes. Clas- 
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Figure 4. Difference in the rest-frame B-band rise time derived 
from the B band, t r {B)(B), and the V band, t r (B)(V), as a 
function of Aniis(B). Differences are generally small, although 
there appears to be a small trend in the rise time as a function of 
Ami5(S). We find evidence for a systematic difference between 
the two measurements. 



sification of SN 1991T/1999aa-like objects requires a spec- 
trum within a week of maximum light to avoid confusion 
with normal SNe la ( ]Li et al.|2001b| ). 

In instances where a SN in our sample has a classifica- 
tion in 



Wang et al. (20091, we adopt their subclassifications. 



Otherwise, our subclassifications are based on BSNIP spec- 
tra. We are left without subclass identifications for only SNe 
1992bo and 1992bc; hence, these two objects are excluded 
from our subsequent analysis. Our final sample consists of 61 
SNe, including 39 normal, 16 HV, and 6 SN 1991T/1999aa- 
like SNe. Figure [2] shows the redshift distribution of our 
sample for each of the three spectroscopic subclassifications. 

We plot the B-band light curves for the 61 SNe that 
passed our cuts in Figure [3] The SNe are shifted relative to 
the time of maximum light and normalised to have the same 
peak flux. Qualitatively, objects with slower post-maximum 
declines (i.e., smaller Am 1 c,(B)) have slower rise times. In 
§4.1, we explore this relationship in more depth. 

H10 analyse a sample of 41 nearby SNe la drawn from 



S07, R99, and the CfA3 sample ( |Hicken et al. 2009a). Our 
analysis benefits from the addition of the LOSS SN la sample 
which has BV data starting at least a week before maximum 
light for 70 SNe. Data from the LOSS sample make up ~70% 
of the objects that pass our cuts. 



3.5 B- and V-Band Results: To Combine or Not 
to Combine? 

In similar analyses of the B-band rise-time distribution, S07 
and H10 combine stretches in the B and V bands using an 
error-weighted mean to produce a final single B-band rise 
time and fall time for each SN. H10 found evidence for a 
weak trend between the difference in stretch values for B 



Figure 5. Difference in the fall times in the B band, tf(B), and 
the V band, tf(V). The difference is consistent with d. 



and V as a function of Arai5(B) using measurements for 
105 SDSS-II SNe la. 

In this paper, we take a different approach and inde- 
pendently measure the B-band and V-band rise times from 
the corresponding photometric data. However, we can use 
the results of our fitting routine to look for similar trends 
in the B-band rise time. In §3.5. 1| we compare the B-band 
rise time derived from the B band to that derived from the 
V band. The B-band rise time found using the V-band data 
is measured by taking the V-band rise time and subtracting 
the time between maximum light in the V and B bands. 

Comparisons by H10 are done in stretch space, while 
ours are between the measured rise/fall time. H10 use a sin- 
gle template rise time of 16.5 d for both the B and V bands 
for fitting the rise stretch of their data light curves. Going 
from rise stretch to rise time only requires multiplying the 
measured rise stretch by the rise time of their fiducial tem- 
plate. Our fitting routine measures the rise and fall values 
of stretch using a Monte Carlo simulation which randomly 
chooses a template rise time based on the best-fit template 
rise time and its Gaussian uncertainty (see j |3.3| for details). 
The measured rise-stretch values are tied to the template 
rise time used for the fit, making the rise time a more ap- 
propriate quantity to compare. 



3.5.1 B-band Rise-Time Comparison 

In Table [l] we show the differences in B-band rise time 
derived using the B band, t r (B)(B), and the V band, 
t r (B)(V), for our sample divided by spectroscopic subclas- 
sification. 

For our entire sample of 61 SNe, we find a mean dif- 
ference of —0.91 ±0.10 d (standard error of the mean), in 
the sense that the B-band rise time is shorter using B-band 
photometry compared to that derived using l/-band pho- 
tometry. Breaking our objects by spectroscopic subclassifi- 
cation, the difference in B-band rise time found using B- 
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Table 1. Mean differences in B-band rise time derived using B- 
and V-band photometry by spectroscopic subclassification 



Subclassification 



t r (B)(B) - t r (B)(V) 



Normal 

High Velocity 

SN 1991T/1999aa-like 

All 



-0.79 ±0.13 d 
-1.08 ±0.16 d 
-1.17 ±0.45 d 
-0.91 ±0.10 d 



Table 2. Mean differences in B- and V-band fall times by spec- 
troscopic subclassification 



Subclassification 



t f (B)-t f (V) 



Normal 

High Velocity 

SN 1991T/1999aa-like 

All 



-0.20 ±0.08 d 
-0.19 ±0.22 d 
-0.32 ±0.39 d 
-0.21 ±0.09 d 



and V-band photometry is —0.79 ± 0.13 d for normal SNe, 
-1.08 ± 0.16 d for HV SNe, and -1.17 ± 0.45 d for SN 
1991T/1999a-like objects. We find evidence for a system- 
atic difference of ~ 1 d between the 2 determinations of the 
B-band rise time. 

In Figure HI we plot difference in B-band rise time as 
a function of Amis(B). There appears to be a small trend 
in the B-band rise time difference with increasing Amis(B). 
Fitting a line results in a slope of — 0.54±0.44 d mag -1 . The 
slope is computed by bootstrap resampling our sample to 
give 1000 realizations of our dataset. The mean and standard 
deviation of the distribution of the fit slopes are adopted as 
the most probable slope value and lcr uncertainty. We do 
not find evidence of a significant trend. 

Our comparison shows evidence for a slight systematic 
trend that the rise time is ~ 1 d longer when measured with 
the V band. At least part of such a shift can be hidden by 
uncertainties in measuring the time of maximum flux of a 
light curve where the derivative slowly approaches within 
a day of maximum light. Combining the two measurements 
will introduce systematic errors into an analysis. 



3.5.2 Fall-Time Comparison 

In a comparison of fall-time stretches derived from the B 
and V bands by H10, the authors found that objects with 
small Ami5(_B) had a larger B-band fall stretch than V- 
band fall stretch (i.e., tf(B) > tf(V)). We look for a similar 
trend in our data. 

Recall that the B-band fall time is defined as the 
amount of time it takes for the B-band light curve to fall 
by 1.1 mag after maximum light in B, and the V-band fall 
time is defined as the amount of time it takes for the V-band 
light curve to fall by 0.66 mag after maximum light in V. By 
construction, our templates both have a fall time of 15 d. 

In Table [2] we show the difference in fall times between 
the B and V bands (i.e., t/(B) — t/(V)). A comparison of 
the fall times for all 61 objects gives a mean difference of 
—0.21 ± 0.09 d (standard error in the mean) between the B 
and V bands for all objects. There is no significant mean 
difference in the fall time measured between the B and V 
bands across all subclassifications. 

We plot the difference in fall times as a function of 
Ami5(B) in Figure [5] The mean fall-time difference is con- 
sistent with d. Fitting a line to the data, we measure a 
slope of —0.47 ± 0.36 d mag -1 , finding no significant evi- 
dence for a trend. 



Table 3. Rest-frame, fall-stretch corrected B- and V-band rise 
times 



Subclassification 


tr'(B) 


tr'(V) 


Normal 

HV 

SN 1991T/1999aa-like 


18.03 ±0.24 d 
16.63 ±0.29 d 
18.05 ± 0.69 d 


20.23 ±0.44 d 
19.43 ±0.33 d 
20.00 ±0.68 d 



4 ANALYSIS 

In this section we present our analysis of the B- and V- 
band rise-time distribution of our nearby sample. Given the 
results found in §3.5| we treat the B and V bands separately. 
In the following analysis, the rise time for a given band is 
defined as the elapsed time between explosion and maximum 
light in that band. Note that the results are not corrected 
for stretch unless indicated by a ' symbol, in which case the 
quantity is corrected for fall stretch. 



4.1 Rise Times Correlated with Decline 

In Figure [6l we plot the rest-frame B-band rise time as a 
function of Amis(B) and the rest-frame V-band rise time as 
a function of Atois(V) for our nearby sample of 61 objects. 
In both bands, there is a strong correlation between the 
rise time and the post-maximum decline, although the scat- 
ter indicates that the situation is more complicated than a 
simple one-to-one mapping between light-curve decline and 
the rise time. In general, SNe with slower post-maximum 
declines (e.g., small Amis(B)/ Amis(V)) have longer rise 
times. This correlation is evident in B, but shows more scat- 
ter in V. Using 61 objects, we calculate a Pearson correlation 
coefficient of -0.69 ±0.03 for B and -0.49 ±0.04 for V. The 
probability of a Pearson coefficient < —0.49 for two uncorre- 
cted variables and 61 measurements is < 0.01%, indicating 
a highly significant correlation in both B and V. The correla- 
tion remains strong even when excluding SN 1991T/1999aa- 
like objects. 

Figure |6] indicates that SN 1991T/SN 1999aa-like SNe 
(i.e., overluminous SNe la) have the longest rise times. In 
contrast, H10 find that overluminous events have some of 
the shortest B-band rise times in their sample. However, 
H10 have few overluminous events and spectroscopic identi- 
fications were not included in their analysis. Instead, iden- 
tifications of overluminous objects were based on Amia(B). 

In Table [3] we give the fall-stretch corrected rise times, 
t r ' , for the various spectroscopic subclassifications in our 
sample for the B and V bands. When B-band light curves 
are fall-stretch corrected to a light curve with Amis(B) = 
1.1 mag (i.e., Sf — 1), the rise time of normal SNe la and HV 
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Figure 6. The rest-frame B-band rise time plotted as a function of h.m\^{B) (left), and the rest-frame V-band rise time as a function 
of Ami^(V) (right). Note that the rise time is not stretch corrected. Using the two-stretch fitting method, we find a correlation between 
Ami5(B) (calculated from sj{B) and EquationHl and rise time. SNe with smaller Ami5(B) (i.e., a slower post-maximum decline rate) 
have longer rise times. SNe that have been identified as HV are plotted as red stars. For fixed Amis(.B), HV SNe appear to have shorter 
rise times than their normal counterparts in the B band. Overluminous SN 1991T/1999aa-like objects are plotted as blue squares. These 
objects have the smallest values of Amis(B) and the longest rise times. A less prominent, but similar correlation exists in Ami5(V) 
(calculated from Sf(V) and Equation |5l and the V-band rise time. 



objects is (respectively) 18.03 ± 0.24 d (uncertainty in the 
mean) and 16.63 ± 0.29 d, a ~ 3<j difference. H10 found an 
average B-band rise of 17.38 ± 0.17 d for the SDSS sample, 
consistent to 2.2<r with our rise time for normal SNe la. 

When we correct our V-band rise times to Amis(V) = 
0.66 mag by dividing by the V-band fall stretch for normal 
objects, we find a fall-stretch corrected, rest-frame V-band 
rise time of 20.23 ± 0.44 d (uncertainty in the mean). The 
fall-stretch corrected rise time of HV objects is 19.43±0.33 d, 
consistent with the V-band rise time of normal objects. If 
we assume that a normal SN la with Amis(_B) = 1.1 mag 
corresponds to Amis(V) = 0.66 mag, then the fall-stretch 
corrected V-band rise time is 2.20 ± 0.50 d longer than 
the B-band rise time. This is within 1.5a of the 1.5 d 
difference reported by R99 despite our disagreement on 
the actual measured B-band and V-band rise times for a 
Amis(B) = 1.1 mag SN (R99 measure 19.5 ± 0.2 d and 
21 ± 0.2 d for their B- and V-band rise times, respectively). 
For the HV objects, we measure a fall-stretch corrected rise- 
time difference of 2.80 ± 0.44 d. The larger difference (in 
absolute terms) in the B- and V-band rise times for HV 
objects compared to normal objects indicates that HV ob- 
jects have a faster rise in the B band than in the V band 
(compared to normal objects). 

4.2 The Rise Time of High Velocity SNe 

Plotted as red stars in FigureRjlare objects spectroscopically 
subclassified as HV, while normal SNe la are plotted as filled 
circles. In the B band, HV objects appear to lie along a locus 
of points below that of normal SNe. For a fixed Ami5(B) 
in the B band, HV objects have a shorter rise than normal 
SNe. Analysing smaller samples, |Pignata et al.| ( |2008[ > and 
Zhang et al. ( 2010 I previously found that HV SNe la appear 



to have a faster B-band rise for a given Ami5(B). This result 
is not evident in the V band, where HV objects do not differ 
significantly from their normal counterparts. 

The ~ 3a difference in fall-stretch corrected rise time 
between HV and normal objects in B and the lack of a 
significant difference in V provides evidence for subtle dif- 
ferences in the photometric evolution of these two subclas- 
sifications. Recent work by |Wang et al.| ( |2009j ) and Foley 
& Kasen (20111 has shown that the two subclassifications 
have different B max — V max pseudocolouJjdistributions, with 
HV objects typically having a redder pseudocolour. Further- 
more, these studies have shown that separating HV and nor- 
mal objects in cosmological analyses reduces the scatter in 
a Hubble diagram from 0.18 mag to 0.12 mag. 



Foley & Kasen (20111 provide a model which offers 



a possible explanation for why HV objects have a redder 
B max — V m ax colour at maximum light. The two dominant 
sources of opacity in the atmosphere of a SN are electron 
scattering and line opacity from Fe-group elements. Elec- 
tron scattering opacity is wavelength independent, while line 
opacity from Fe-group elements is most significant at wave- 
lengths shorter than 4300 A ( |Kasen fc Plewa|2007[ hereafter 
KP07). The transition from electron scattering to Fe-group 
line opacity occurs near the peak of the B band. SNe with 
high-velocity ejecta will have broader Fe-group absorption 
features, which will decrease flux in the B band while having 
little effect on V-band flux. 



3 Note that B max — Vmax is the ZJ-band magnitude at maximum 
light in the B band minus the V-band magnitude at maximum 
light in the V band. The dates of maximum light are typically 
offset by ~2 d, making this quantity not an actual observed colour 
of the SN at a discrete time. 
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Foley & Kasen (20111 analyse models of an off-centre 



failed deflagration to detonation from KP07 to explore the 
expected differences in observables of HV and normal SNe. 
The KP07 models studied a single SN with an off-centre 
ignition viewed from different angles. When viewed on the 
side nearest the ignition, the KP07 models produce a SN 
with HV features. When viewed from the side opposite the 
ignition, we will view a normal SN. Although the models 
presented in KP07 were intended to study the observational 
consequences of off-centre explosions, their model spectra 
and light curves conveniently produce a similar distribution 
of ejecta velocities to that observed between normal and HV 
SNe depending on viewing angle. The set of models predicts 
that HV objects should have redder colours at maximum 



light in comparison to normal objects (see Fig. 8 of Foley & 
KasenpOll \ on theoretical grounds. 



Using the model light curves of KP07, we can explore 
how the rise time of HV objects will differ from that of 
normal objects. If all other parameters are equal (i.e., Nickel 
mass, kinetic energy, etc.), the KP07 model predicts that HV 
SNe should have shorter rise times in B compared to normal 
SNe (as indicated by the model light curves in their Fig. 11) 
due to the enhanced line opacity from Fe-group elements. 
Although not specifically stated by KP07, one would expect 
the V-band rise time of the two subclasses to be similar, 
since the opacity in this wavelength range is dominated by 
electron scattering. 

Consequently, the KP07 models also predict that the 
enhanced opacity at blue wavelengths from HV components 
in the SN ejecta reduce the peak absolute _B-band magni- 
tude and hasten the evolution of the _B-band light curve, 
producing a larger Amis(B). This complicates a direct com- 
parison between the rise time of HV and normal SNe and 
does not necessarily predict the observed result in the left 
panel of FigurelS] that HV SNe lie on a locus of points below 
that of normal SNe in the Ami$(B)-t r plane. One expec- 
tation of the KP07 models is a Ami5(B) distribution for 
HV objects that is pushed to larger values in comparison 
to normal objects. We find a median Ami$(B)= 1.09 mag 
for the HV objects and a median Amis(B)— 1.11 mag for 
the normal objects. However, our sample is most likely not 
representative of a complete sample of SNe la and suffers 
from observational bias. 

Without quantifying how HV ejecta change both 
Ami5(_B) and t r , we cannot definitively claim that the mod- 
els of KP07 explain our result. While it is beyond the scope 
of this paper to match theoretical models to our observa- 
tions, the KP07 models offer a possible qualitative expla- 
nation for why there are significant differences between HV 
and normal SNe in the B band, but not in the V band. 



Recent work by Maeda et al. (20101 finds evidence that 



SNe with a high velocity gradient (HVG) in the Si II line 
may be the natural result of viewing an asymmetric ex- 



plosion. HVG SNe are defined by Benetti et al. ( 2005 1 as 



SNe that exhibit a time derivative in the velocit y of the 
1 d _1 around maximum light. Benetti 



Si II line > 70 km s _i d _1 around maximum light 
|et al.| provide evidence that the HV and HVG subclassifi- 
cations are highly correlated. Looking at the Fe II A7155 
and Ni II A7378 nebular emission in late-time spectra that 
trace the deflagration ashes, Maed aet al.| find that HVG SNe 
tend to exhibit redshifted lines while low-velocity gradient 
(LVG) SNe show a blueshift. The authors attribute this ob- 



servational distinction to the difference between viewing an 
asymmetric explosion from the side nearest the site of initial 
deflagration (LVG) and the opposite side (HVG). The mod- 
els of asymmetric explosions of SNe la presented by |Maeda| 
et al. (2011 hereafter Mil) indicate that viewing a SN from 
the far side, one expects HVG SNe to have longer bolomet- 
ric rise times and smaller Amn,(B) than comparable LVG 
SNe. This effect is more prominent in SNe with less 56 Ni 
(for instance, compare models A0.3 and A0.6 in Fig. 11 of 
Mil). The longer rise times in HVG objects is attributed to 
an increased optical depth for optical photons when viewing 
the explosion from the side opposite the explosion site (i.e., 
the site of 56 Ni synthesis) . 

This is not necessarily inconsistent with our observa- 
tional result that the HV SNe have a different rise-time dis- 
tribution than normal SNe. If we view the SN from the side 
nearest the initial explosion site, we will see a LVG SN, 
and if we view the SN from the opposite side, we will see 
a HVG SN. Assuming HVG implies HV, we will measure a 
longer rise time and a smaller Ami$(B) for a high- velocity 
object assuming the B band roughly traces the bolomet- 
ric behaviour. This will move the HV object up and to the 
left of a normal SN in the Ami*,(B)-t r plane (Fig.pl. De- 
pending on how the viewing angle affects both rise time and 
Ami5(B), the models of Mil could put the locus of HV 
points below that of normal SNe. However, the models pre- 
sented in Mil do not outright predict a difference in rise 
times of HV objects in the B and V bands. 

The models of KP07 and Mil offer opposite theoret- 
ical predictions for the rise time of HV SNe in compari- 
son to normal SNe. KP07 predict that HV SNe should have 
larger Amis(B) and shorter rise times than normal SNe, 
while Mil predict smaller Am\$(B) and longer rise times. 
The differences are rooted in the nature of the SN asym- 
metry (i.e., the distribution of intermediate-mass elements 
and 56 Ni) and the treatment of opacity. KP07 use an ex- 
pansion opacity formalism which sums over individual lines, 
while Mil use a frequency averaged gray opacity. A better 
test of the HV models may possibly be found by looking 
at the Amis(B) distributions of a complete SNe la sample, 
since both models predict that that asymmetries should in- 
fluence the measured Am^(B). |Wang et al.| ( |2009[ ) found 
the Amis(B) distribution of HV and normal SNe to be 
strikingly similar despite different _B max — Kn ax distributions, 
although it is not clear that their sample is complete. 

Despite evidence that HV objects have different rise- 
time properties than normal SNe, a clear physical picture 
remains elusive. Further efforts in modeling HV objects may 
shed light on the rise-time distribution of the different spec- 
troscopic subclassifications. 

4.3 The Rise Minus Fall Distribution 

In the top two panels of Figure [7] we compare the rise time 
minus the fall time (RMF) as a function of Amis(B) for 
the B band (left panel) and Amis(V) for the V band (right 
panel). Note that the RMF has not been stretch corrected. 
As in Figure [6] blue squares refer to SN 1991T/SN 1999aa- 
like objects, red stars refer to HV SNe la, and black circles 
refer to spectroscopically normal SNe la. Overplotted in a 
broken line is the expectation from a one-to-one mapping of 
rise time to fall time using a single-stretch parametrisation. 
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Figure 7. Top panels: The rest-frame corrected rise time minus the fall time (RMF) as a function of h.rn\§(B) for the B band (left) and 
of Ami5(V) for the V band (right). These quantities are not stretch corrected. Plotted as a dashed black line is the prediction using a 
single-stretch prescription where there is a one-to-one mapping of rise time to fall time. The nonlincarity in the prediction is due to the 
quadratic relationship between Sf and Amis(B) in EquationRI There is significant scatter about the line. We find that slowly declining 
SNe tend to have a faster rise than predicted by the single-stretch model. Bottom panels: The RMF distribution in the B and V bands. 
HV objects have a shorter B-band RMF in comparison to objects which are spectroscopically normal. A K-S test indicates a < 0.01% 
chance that the two distributions are drawn from the same parent distribution. No significant difference is seen in the V-band RMF 
distribution. SN 1999ee is the only normal SN with B-band RMF < — 1 d, although a previous study noted HV features in spectra of 
this object (|Mazzali et al.|2005|). 



For the B band we use a fiducial 18.03 d rise time and for 
the V band we use a rise time of 20.23 d based on the fall- 
corrected rise times found in §4.1| Clearly, our sample does 
not strictly follow a single-stretch parametrisation. Similar 
to the results of H10, a number of slowly declining objects 
have a faster rise time (i.e., a smaller RMF) than expected 
from a one-stretch parametrisation in both B and V. We re- 
iterate that based on the results of Figure [6] more luminous 
SNe have longer rise times than less luminous SNe; however, 
Figure [7] indicates that more luminous SNe have faster rise 
times than expected based on a single-stretch parametrisa- 
tion. 

In the bottom two panels of Figure [7] we plot the rest- 
frame B- and V-band distributions of RMF for the various 



subclass identifications. In B, HV objects have a mean RMF 
of 1.55 ± 0.27 d (uncertainty in the mean) in comparison 
to normal SNe that have an RMF of 2.77 ± 0.20 d. The 
4<t difference in mean RMF offers evidence that the two 
distributions may be drawn from different populations. The 
mean RMF for our sample by subclassification can be found 
in Table p] Despite significant differences in the B band, the 
mean V-band RMFs for the spectroscopic subclassifications 
are consistent with one another. 

A Kolmogorov-Smirnoff (K-S) test finds that the two 
groups of SNe have a ~ 0.01% probability of being drawn 
from the same B-band RMF distribution and a ~ 44% prob- 
ability of being draw from the same V-band RMF distribu- 
tion. The lone SN with RMF < -1 d in B is SN 1999ee. 
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Table 4. Mean rise minus fall times for different spectroscopic 
subclassifications 



Subclassification 




RMF (B) 


RMF(V) 


Normal 
High Velocity 
SN 1991T/1999aa 


-like 


2.77 ±0.20 d 
1.55 ±0.27 d 
3.38 ±0.73 d 


4.83 ±0.33 d 
4.45 ± 0.28 d 
5.72 ±0.73 d 



Note - Quantities have not been corrected for stretch. 



Mazzali et al. ( 2005 1 find evidence of high- velocity features 
in spectra taken before maximum light. However, we retain 



the classification of Wang et al. ( 2009 1 and regard SN 1999ee 
as being normal. 

Our application of the K-S test indicates that HV and 
normal objects have a high probability of being drawn from 
different populations. However, the K-S test does not reveal 
whether this difference is physical in origin or a reflection of 
observational bias. For instance, the difference in RMF be- 
tween the two populations may also be a result of different 
stretch distributions. Such differences cannot be disentan- 
gled without knowledge of the observational bias in each of 
the photometric surveys used in our sample. Future surveys 
with proper spectroscopic follow-up observations and under- 
standing of biases will aid in exploring the difference in RMF 
distribution for HV and normal SNe. 

Our sample also includes 6 objects which show spec- 
troscopic similarities to the overluminous SN 1991T or SN 
1999aa. Based on the B- and V-band distribution of RMF 
for SN 1991T/1999a-like objects, these objects have slightly 
larger RMFs compared to normal and HV objects, however, 
the mean of the 1991T/1999aa-like distribution is consistent 
with the normal objects within lcr. 

S07 find evidence for two distinct RMF distributions 
using eight SNe. While we have presented evidence for dis- 
tinct HV and normal populations, S07 only have a single 
HV object in their sample (SN 2002bo). The other seven 
SNe are classified as normal by the criterion of Wan g et al.| 
( 2009 1 . We do not see evidence for two populations of rise 



times within our spectroscopically normal SNe. 

Similar to the results of H10, we find that our sam- 
ple does not follow the expectation from a one-stretch 
parametrisation of light-curve shape. This is especially ev- 
ident in V, where the trend appears to go in the opposite 
direction predicted by a single-stretch fit. 



4.4 Rise-Time Power Law 

Previously, we assumed the rise in flux took the form of a 
parabola, n = 2, based on physical arguments presented in 
§3.1. In this section, we fit for the functional form of the 
rising portion of the light curve (r ^ —10 d) as a power 
law of the form / = A(r + t r ) n , where t r is the rise time. 
Allowing n to vary, we perform a \ 2 minimization to find 
the best-fit power law to our B-band photometry. However, 
unlike the analysis in previous sections, we restrict our two- 
stretch fitting procedure to — 10 < r < +35 d in order to 
avoid imposing a shape on the region we plan to fit. We find 
a best fit of n — 2.20l 019 , consistent (lcr) with the expand- 
ing fireball modeled by a parabola. The uncertainty in the 
power-law index is found using a Monte Carlo simulation 



similar to that outlined in §373] except modified to only fit 
within the region — 10 < r < 35 d. 

Our result is in agreement with that of C06, who find 
n = 1.8 ± 0.2, and H10, who find n = 1.80± O if. However, 
H10 show evidence of significant colour evolution during this 
period of light-curve evolution, challenging our assumption 
of modest temperature change. They find a linear B — V 
colour evolution of 0.5 mag between 15 d and 9 d before 
maximum light, leading to a temporal dependence closer to 
n ~ 4 rather than the n — 2 predicted by an expanding 
fireball. 

We check for similar evidence of colour evolution in the 
expanding fireball phase in our sample. Using the colour 
curves for SNe in the range of —15 < r < —9 d relative 
to maximum light, we create a median B — V colour curve 
similar to the analysis presented by H10. We divide the data 
into 1 d interval bins and take the median B — V colour for 
each bin. Corrections for Milky Way extinction were made 
using the dust maps provided by Schlegel et al. (19981. No 



corrections were attempted for possible host-galaxy extinc- 
tion. Our median colour curve shows a small change of 0.05 
mag over the 6 d interval. However, inspecting the B — V 
colour curve for SN 2009ig in this range, we find a drop 
of 0.5 mag over the same interval, similar to what is re- 
ported by H10. This leads us to suspect that our median 
colour curve does not necessarily reflect the full sample, and 
at least in individual SNe, there can be significant colour 
change at r < —9 d. The discrepancy between the median 
colour curve and that of SN 2009ig may be a result of using 
noisy data in the earliest time bins for the median colour 
curve. 



4.5 Companion Interaction 



Recently, Kasen (20101 proposed using emission produced 



from the collision of the SN ejecta with its companion star 
to probe the progenitor system. Similar to the shock break- 
out in core-collapse SNe that is theoretically well understood 



(e.g. , |Klein fc Chevalier|1978~l|Matzner fc McKee|1999[) and 
observed (e.g., |Soderberg et al.|2008||Modjaz et al.|2009| ), in 
the single-degenerate progenitor scenario the expanding SN 
ejecta are expected to collide with the extended envelope of 
the mass donating companion. In the case of a companion 
undergoing Roche-lobe overflow, the radius of the compan- 
ion is on the same order of magnitude as the separation 
distance. The timing and luminosity emitted from the inter- 
action will depend on the mass of the companion and the 
separation distance. The shock emission is expected to be 
brightest in the ultraviolet, but detectable in the B band. 



Hayden et al. ( 2010b I looked for this signal in the rising por- 



tion of the B-band light curves of 108 SDSS SNe la, finding 
no strong evidence of a shock signature in the data. Using 



simulated light curves produced with SNANA (Kessler et al. 



2009) and a Gaussian to model the shock interaction, the 
authors constrain the companion in the single-degenerate 
scenario to be less than a 6 M© main-sequence star, strongly 
disfavouring red giants (RGs) undergoing Roche-lobe over- 
flow. 

For this analysis, we focus on our B-band data where 
we have the best chance of detecting the signs of shock emis- 
sion. |Kasen| p010| predicts that in the initial few days af- 
ter explosion, the luminosity produced by the interaction 
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Figure 8. Two-stretch corrected B-band light curves compared to the expected flux of a collision between the SN ejecta and a 1 Mq 
companion assuming the companion is undergoing Roche-lobe overflow. Only data in the range —10 < r < +35 d are used to stretch 
correct the light curves, to avoid imposing a shape onto data at r ^ — 10 d. Plotted as a red line is the light-curve template with no 
shock interaction assuming various power laws for the initial rise, indicated by n, and rise times, t r . Plotted as a blue dashed line is 
the expected "shocked template," which is the expected shock emission added to the template assuming different separation distances 
013 and the fraction of SN explosions that show signs of interaction, /. The expected flux from the companion is estimated using the 
analytic models of |Kasen| l |2010[ | . The bottom plot for each panel shows the resulting normalised residual curve between the data and 
the "shocked template." We find a considerable degree of degeneracy between the adopted rise time, the power-law index for the initial 
rise, and the shock emission. As the four plots indicate, a similar minimum reduced x 2 can be achieved with n = 2, t r = 17.92 d and no 
shock emission (top left) or by varying the rise time, power-law index, and amount of shock emission (top right, bottom left, and bottom 
right, respectively). 



with the companion will dominate the luminosity powered 
by 56 Ni decay. Inspection of individual light-curve fits with 
the two-stretch fitting routine do not show the tell-tale signs 
of strong interaction. To take advantage of the power in 
numbers, we analyse all of the light curves as an ensemble. 
We then compare data from the earliest light-curve epochs 



to the models of Kasen (2010) to place constraints on the 



mass and distance to the companion. 

We start by applying our two-stretch fitting routine 



again, but limiting the fit to data — 10 < r < +35 d in 
order to avoid forcing a shape on the earliest data and pos- 
sibly suppressing the signature of interaction. We construct 
models of "shocked" template light curves including contri- 
butions from companion interaction using the analytic so- 
lutions for the properties of the emission found by [Kasen] 
(|2010[) . Using equations for L Cl j so and T e g from Kasen (2010 1 



(Eq. 22 and 25, respectively) and assuming the emission is 
that of a blackbody, we calculate the expected flux density 
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at the peak of the B band (A e fi = 4450 A) at a distance of 
10 pc. We normalise the flux density from the shock to a 
peak SN la magnitude of Mb = —19.3 mag. 

Based on the opening angle of the shock interaction, 



Kasen (20101 predicts that the interaction signature should 
be visible in ~ 10% of SN la explosions. Ideally, one should 
account for the viewing angle dependence of the detected 
shock emission. We make the simplification that we are 
either looking along the axis or we are not. Our light- 
curve models including the effects of interaction are the un- 
shocked template plus a fraction, /, of the collision flux cal- 



culated using the analytic model of Kasen (20101. A model 
with / = represents a situation in which shock emission 
is never detectable (and reduces to the unshocked template) 
and / = 1 is a scenario in which shock emission is detectable 
in every SN explosion. 

For most of the analysis in this paper, we adopted an 
expanding fireball model for the early rise of the light curve 
(r ^ —10 d) which assumed a power-law index of n = 2. 
Under this assumption, we found that our early-time data 
were best fit by a rise time of 17.92 d. However, as discussed 
in §4.4[ the assumption that the initial rise of the light curve 
has n = 2 may be somewhat questionable given rather sig- 
nificant changes in SN colour at early phases. As a test of the 
degeneracies between the added shock emission and the as- 
sumptions that go into constructing an unshocked template, 
we consider a number of different templates with different 
rise times and power-law indices. In addition to our nominal 
template with t r — 17.92 d and n = 2, we also try templates 
with t r — 17.11 d and n = 2, t r = 17.92 d and n = 2.2, and 
t r = 17.11 d and n - 1.8. 

The parameters free to vary in the model for shock in- 
teraction are 0,13, the distance to the companion normalised 
to 10 13 cm, M c , the companion mass, and /, the fraction of 
SN explosions that produce detectable shock emission. For 
this simple analysis the probed mass is fixed at M c = 1 Mq, 
to explore RGs as a possible companion, and we set / equal 
to either (i.e., no emission is detected) or 0.1, the expected 
fraction of SNe with detectable shock emission based on the 
opening angle of the shock. For each of our four different 
unshocked templates, we calculate the minimum \ 2 statis- 
tic for models by varying a.13 to find the best-fit shocked 
template. The fit is restricted to data within the range of 
-16 ^tsC -10 d. 

We note that the minimum reduced x 2 in all of our fits 
for ai3 exceeds 1 and is usually closer to 3-4. This is not 
completely unexpected since data at r ^ — 10 d were not in- 
cluded in the two-stretch fitting procedure used to normalise 
the light curves, and this may induce correlated errors into 
the data at r ^ — 10 d. We also suspect that the reported 
errors for the earliest photometry epochs may be underesti- 
mated. Furthermore, the scatter in the data points may even 
be the physical result of a distribution of companion masses, 
separation distances, and viewing angles contributing to the 
emitted flux from shock interaction. Unfortunately, we are 
not in a position to disentangle what is contributing to the 
scatter. 

Our final results for each template are shown in Figure 
[8] Plotted in red are the unshocked templates and in blue 
are the shocked templates including some contribution of 
companion interaction. In each case of an unshocked tem- 
plate, we can find an acceptable fit with similar \ 2 by vary- 



ing the separation distance. There is a significant degree of 
degeneracy between the parameters, making the task of dis- 
entangling the true sign of companion interaction at this 
level extremely difficult. For instance, for our template with 
n — 2 and t r — 17.92 d, we do not see signs of any inter- 
action. However, had we used a template with n — 2 and 
t r — 17.11 d, we would find that if the fraction of SNe that 
showed signs of shock interaction is / = 0.1, then % 2 ls min- 
imised with ai3 = 0.9. Similarly, when varying n and t r in 
the unshocked template, suitable matches to the data can 
be found by adjusting the amount of shock interaction. 

The degeneracy is partially broken by increasing the 
contribution of companion interaction to the "shocked tem- 
plate" by either increasing / or 013. Fixing / = 0.1 based on 
the expected opening angle of the shock, a.13 = 2 increases 
the reduced x 2 by more than 1 for all of our unshocked 
templates. 

In summary, we find that with our limited sample of 
early-time data, we cannot completely disentangle the be- 
haviour of the power-law index, the adopted rise time, and 
the amount of companion interaction. This result makes use 
of all of the light-curve data at r ^J — 10 d as an ensemble. 
Better constraints could be placed on individual, excellently 
sampled light curves in the rise region. In particular, efforts 
to obtain data at bluer wavelengths, where the shock inter- 
action is stronger and more easily detectable, will greatly 
help to break degeneracies in the fitting process. 



Recently, Justham (20111 and Di Stefano et al. (20111 



have argued that the transfer of angular momentum from 
the companion donor star to the C/O white dwarf acts to 
increase the critical mass required to explode and the time 
required for the white dwarf to explode. Consequently, the 
donor star may possibly evolve past the red giant phase, 
vastly reducing the cross section for interaction with SN 
ejecta. The subsequent shock emission for a RG would de- 
crease beyond current detection limits, thus saving RGs as 
a possible donor in the single-degenerate scenario. The in- 
creased time to explosion allows the mass ejected from the 
envelope of the donor star time to diffuse to the density of 
the ambient interstellar medium, offering an explanation for 
the tight constraints placed on the presence of Ha in nebular 
spectra of SNe la ( |Leonard|2007| and references therein) . 

4.6 SN 1991bg-Like Objects 

The analysis presented thus far has excluded SN 1991bg- 
like objects due to fits which had unacceptably large x 2 P er 
degree of freedom. To investigate what was causing the poor 
fits, we focus on the B-band light curve of the SN 1991bg-like 
SN 1999by from the LOSS sample, which has data starting 
about — 10 d before maximum light in the B band. As seen in 
Figure [9] the largest difference in light-curve shape between 
SN 1999by and our template occurs at t > +10 d, where 
the light curve of SN 1999by transitions to a slower linear 
decline not seen until t > +30 d in normal SNe la. Similar 
differences in light-curve shape for SN 1991bg are found by 



Filippenko et al. (1992b I and Leibundgut et al. (19931. The 



fit parameters of s/ and s r are correlated; thus, a bad fit 
for Sf will propagate into an incorrect determination of s r . 
Given the different light-curve shapes for SN 1991bg-like 
objects and normal objects for t > 10 d, we are unable to 
acceptably fit Sf and thus s r . For the purposes of exploring 
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Figure 9. Comparison of the B-band light curve for the SN 
1991bg-like SN 1999by (solid circles) and the LOSS template 
(black dashed line) stretched along the time axis by 0.74 to match 
the rise portion of the light curve. The light-curve evolution of SN 
1999by for t > 10 d past maximum makes it impossible to find a 
reasonable fit using our two-stretch fitting routine. If we restrict 
our fit to the pre-maximum portion of the light curve, we find a 
rise time of 13.33 ± 0.40 d. 



the rise time of SN 1991bg-like objects, we restrict our fit 
to the pre-maximum portion of the light curve using the 
date of maximum found with a low-order polynomial fit. 
For SN 1999by, we find a best-fit rise time of 13.33 ± 0.40 d 
(overplotted as dashed lines in Fig. [9| , indicating that they 
join other objects with large Amis(B) as the fastest risers 
in our sample. This matches the qualitative results found by 



Modjaz et al. (20011 for the SN 1991bg-like SN 1998de. 



4.7 Impacts of Fitting Cuts 

To ensure that the results described in the previous sections 
are not a result of the cuts made in §3.4| we reanalyse our 
data, both tightening and relaxing the constraints on the 
uncertainty in t r and tj , <Jt r t ■ Most of the reported results 
are not highly sensitive to fitting cuts. Restricting accept- 
able fits to reduced \ 2 = 1-5, a first epoch at r = —10 d, 
and at r . — 1 d decreases the number of available ob- 
jects to 26 normal, 9 HV, and 3 SN 1991T/SN 1999aa-like 
SNe la. The probability that the B-band RMF populations 
are drawn from the same parent population increases to 
~ 4%. The difference in the B-band fall-time corrected rise 
time, t r ', between HV and normal SNe remains significant at 
— 1.44±0.49 d, indicating a faster rise for HV objects. Relax- 
ing constraints to a reduced % 2 — 2, first epoch at r = —5 d 
relative to maximum light, and at,, f — 1-5 d increases the 
difference in fall-stretch corrected B-band rise time between 
HV and normal objects to —1.86 ± 0.49 d. The U-band rise 
times remain consistent between HV and normal objects, as 



do the U-band RMF distributions. Overall, changing what 
we define as an "acceptable" fit does not impact our results. 



5 DISCUSSION 

We have presented an analysis of the rise-time distribution 
of nearby SNe la in the B and V bands. Using a two- 
stretch fitting technique, we find that the SN rise time is 
correlated with the decline rate in the sense that SNe with 
broader light curves post-maximum (i.e., light curves with 
small Amis (-B)/ Amis (V)) have longer rise times. While SN 
1991bg-like objects could not be fit well by our two-stretch 
fitting procedure, we found that restricting our analysis to 
the pre-maximum data for the SN 1991bg-like SN 1999by 
gives a rise time of 13.33 ± 0.40 d. This follows the expected 
trend of a fast rise leading to a fast decline. 

Using a sample of 105 SDSS SNe at intermediate red- 
shifts (0.037 < z sC 0.230), H10 find that there is a great 
diversity in B-band rise times for a fixed Amis(B), and 
that the slowest declining SNe tend to have the fastest rise 
times. While we do see evidence for scatter in rise times with 
Ami5(_B), we find a strong correlation between Amis(B) 
and rise time in the sense that slowly declining SNe la have 
longer rise times. This correlation is strong in the B band, 
but weaker in the V band. The discrepancy with H10 may 
be a result of combining B- and U-band stretches which was 
avoided in this analysis. 

We find that a single value of the stretch does not ad- 
equately describe the rising and falling portions of our light 
curves. Similar to the results in H10, more luminous SNe 
with longer fall times have shorter rise times than one would 
expect from a single-stretch prescription for light curves. 
This is especially evident in our analysis of U-band light 
curves. However, we reiterate that while luminous SNe have 
shorter rise times than expected from a single-stretch pre- 
scription, they still have longer rise times than less luminous 
SNe, contrary to the findings of H10. 

R99 and C06 found a fiducial B-band rise time of 
~ 19.5 d for a "typical" SN la. In a departure from previous 
results, H10 report an average rise time of 17.38 ± 0.18 d. 
We find a fall-stretch corrected (i.e., corrected to have a 
post-maximum fall of Amis(B) = 1.1 mag) B-band rise 
time of 18.03 ± 0.24 d for spectroscopically normal SNe and 
16.63 ± 0.30 d for HV SNe. Our B-band rise time for spec- 
troscopically normal SNe with Amis(B) = 1.1 mag is in 
agreement with the average rise time found by H10 at the 
2.2a level. When correcting our U-band light curves to a 
post-maximum fall of Amis(U) = 0.66 mag, we find a fall- 
stretch corrected U-band rise time of 20.23 ± 0.44 d. 

After correcting for post-maximum decline rate, HV 
SNe la have faster rise times than normal SNe la in the 
B band, but similar rise times in U. We find a ~ 3<r dif- 
ference in the fall-stretch corrected rise time between HV 
and normal SNe in the B band. The rise minus fall (RMF) 
distributions (not corrected for stretch) of the two subclassi- 
fications show significant differences in the two populations. 
The peak values of the distributions are offset by 1.22±0.34 d 
in the B band, with HV objects having a faster rise time. 
A K-S test indicates a ~ 0.01% probability that HV SNe 
come from the same parent RMF population as normal SNe 
in the B band. Despite differences in the B band, we see no 



© 2011 RAS, MNRAS 000,[T1|T7] 



16 Ganeshalingam, Li, & Filippenko 



evidence of a difference in the RMF populations in the V 
band. 

Based on the model presented by|Foley fe Kasen (2011 1 



and the models of KP07, we offer a possible qualitative ex- 
planation for why HV SNe should have a different B-band 
rise-time distribution than normal SNe, but similar V-band 
distributions. The physical origin of the difference is possibly 
rooted in the different opacity mechanisms at work in the 
B and V bands. Line blanketing from Fe-group elements is 
the dominant source of opacity for wavelengths shorter than 
~ 4300 A, the peak of the B band. At longer wavelengths, 
such as the V band, electron scattering is the dominant 
source of opacity. Rapidly moving ejecta, as is the case with 
HV objects, will broaden absorption features, diminishing 
the B-band flux without affecting the l/-band flux. Mod- 
els from KP07 show that all other things being equal (e.g., 
Ni mass, kinetic energy), this leads to faster light curves 
with larger Amin(B) for HV objects. However, the enhanced 
opacity at short wavelengths also affects Amis(B), compli- 
cating the application of the models to our result. Further 
modeling and observations of HV SNe fa are required to 
shed light on the photometric differences of this spectro- 
scopic subclass. 

We fit the earliest data in our sample (r ^ —10 d) to 
find that the flux rises as a power law with index 2.20lo'jg. 
This is consistent with the unimpeded, free expansion of 
the expanding fireball toy model that predicts an index of 
2. However, a preliminary analysis of SN 2009ig in the range 
— 15 < r < —9 d shows evidence of significant colour evolu- 
tion, contrary to the assumption of little to no colour evolu- 
tion in the expanding fireball model. H10 find similar colour 
evolution in an analysis of the B — V colour curve of SDSS 
SNe la and derive an expected power-lax index of 4. 

We compare our early-time B-band data as an ensem- 
ble to models (Kasen 20101 of shock interaction produced 



from SN ejecta colliding with the mass-donating companion 
in the single-degenerate progenitor scenario. When relaxing 
our assumptions on the functional form of the early-time 
template light-curve behaviour (i.e., changing the rise time 
or power-law index of the rise), we find that our data re- 
quire some amount of shock interaction to remove system- 
atic trends. This indicates a level of degeneracy between the 
adopted template rise time, the power-law index, and the 
amount of shock interaction required to match the data. 

Future surveys with high-cadence search strategies will 
provide well-sampled SN light curves starting days after ex- 
plosion, substantially adding to the sample of rise-time mea- 
surements. Complemented with spectroscopic follow-up ob- 
servations, analysis of the RMF distribution can be further 
broken into different spectroscopic classes to provide insights 
into the underlying populations and the physics that differ- 
entiates them. 
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